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SECTION  I 
INTRODUCTION 

The  effects  of  external  store  carnage  on  aircraft  performance  can  be 
significant,  especially  for  fighter  aircraft  carrying  fuel  tanks  or  air  to 
ground  weapons.  Degradation  in  top  speed  and  range  must  be  accurately 
estimated  to  determine  if  a  particular  store  configuration  is  practical  from 
a  mission  standpoint.  Flight  testing  can  provide  answers,  but  it  is  extremely 
expensive  and  time  consuming.  For  this  reason,  only  selected  aircraft  loadings 
are  usually  tested. 

To  estimate  performance  parameters  for  a  greater  number  of  configurations, 
wind  tunnel  test  data  are  often  used.  A  scale  model  of  the  clean  aircraft  and  one 
with  the  aircraft  loaded  with  the  store  configuration  of  interest  is  tested. 

Curves  of  lift  versus  drag  coefficient  are  constructed,  and  the  drag  increment 
due  to  external  stores  is  calculated  for  a  particular  flight  condition,  i.e.  , 
lift  coefficient. 

Because  of  the  nature  of  any  measurement  process,  an  uncertainty  is 
associated  with  the  wind  tunnel  data.  The  total  data  uncertainty  is  a  com¬ 
bination  of  uncertainties  in  tunnel  test  conditions,  model  positioning,  and 
model  instrumentation.  Since  the  lift  coefficient  value  for  the  flight  con¬ 
dition  of  interest  rarely  appears  as  a  test  point,  a  fitted  curve  must  be 
constructed  from  the  data.  This  curve  fit  process  is  another  source  of  possible 
uncertainty.  This  report  investigates  the  uncertainty  associated  with  wind 
tunnel  drag  data  and  presents  a  method  that  can  be  used  to  calculate  it  for  a 
specific  flight  condition  of  interest. 
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SECTION  II 


GENERAL  METHOD  OF  UNCERTAINTY  CALCULATION 

Uncertainty  associated  with  C^  at  a  given  condition  is  caused  by 
uncertainties  in  the  measurement  process  and  those  induced  by  the  curve  fit 
procedure.  This  requires  that  the  uncertainties  associated  with  each  factor 
be  propagated  to  arrive  at  a  final  result.  The  method  chosen  here  is  the  Taylor 
series  method  of  error  propagation. 

A  derivation  of  the  Taylor  series  method  can  be  found  in  Reference  1.  Some 
of  the  assumptions  used  in  the  derivation  are: 

1.  Response,  Z,  is  defined  as  a  function  of  the  measured  variables 

Xi.  X2,  ...»  Xn» 

2.  Z  is  continuous  in  the  neighborhood  of  u,  ,  y,  ,  ...,  p,  .  yv  , 

Xj  X2  Xn  Xj 

y  ,  ...,  y  are  the  mean  values  associated  with  Xj,  x2,  ...,  x  which  all  have 
x2  xn  n 

error  distributions  about  the  point  of  interest. 

3.  Z  has  continuous  partial  derivatives  in  the  vicinity  of  u  ,  u  ,  .... 

xi  x2 

\ ' 

4.  Xj,  x2*  ...»  x  are  independent  of  each  other. 

5.  (pv  -  Xj),  (y  -  x2 ) ,  ...,  (y  -  x  )  are  small  or 

X1  x2  11 

32z  ,  a2z  ,  ...,  32z  are  small  or  zero. 

?  2“  2 

3Xi  ax2  ax 

i  ^  n 

The  assumptions  will  be  satisfied  if  the  functions  considered  are  restricted 
to  smooth  curves  near  the  point  of  interest  with  no  discontinuities  and  with  higher 
order  derivatives  either  small  or  zero. 

The  results  of  the  derivation  show  that: 


2 


If 

Z  =  f  (xj ,  x2,  . . . ,  x^)  ( 1) 

then 


S(Z)  ='  {  f-r§-  s<xi)  1  "  +t 


D  X 


3  Z 

3  X, 


S(x2)  ]  2  +  ...+[  ||-  S(xn)  ]  2  } 


(2) 


where  S(Z),  S( X2 ) ,  S(X2),  ...,  S(Xn)  are  the  precision  indices  of  the  response, 
Z,  and  the  variables  Xj,  X2,  X  .  The  precision  index  is  the  computed 
standard  deviation  of  the  measurements  (i.e.,  random  error).  It  is  defined  as 


S  = 


n^T  l  -  *>2 


(3) 


Xi  and  x  are,  of  course,  the  value  of  x  of  a  particular  point  and  the  mean 
value  of  x,  respectively.  If  the  sample  size  is  large,  the  precision  index  is 
approximately  equal  to  the  actual  population  standard  deviation  (a)  associated 
with  the  random  variable  Z. 

Depending  on  the  confidence  level  attached  to  the  response,  the  uncer¬ 
tainty  is  simply  a  function  of  the  precision  index  or  standard  deviation. 

This  assumes  that  no  bias  or  systematic  errors  are  present.  For  drag  increments 
caused  by  external  stores,  this  is  a  reasonable  assumption.  The  drag  increments 
are  calculated  as  differences  between  data  taken  during  one  test  with  the  same 
model  and  instrumentation.  While  bias  errors  may  be  present,  they  are  approximately 

equal  for  different  model  configurations.  When  the  increment  is  determined,  the 
bias  errors  should  drop  out.  This  assumption  is  of  critical  importance  in  the 
following  discussions.  It  would  not  be  true  if  raw  data  from  separate  tests 
were  compared  or  if  different  instrumentation  was  used  during  one  test. 


For  the  case  of  interest,  Cp  and  Cp  are  both  random  variables  since  both 
are  measured  during  testing.  With  the  assumption  that  measured  values  of  Cl 


associated  with  Cp  on  the  fitted  Cp  -  Cp  curve  is  found  as  follows. 
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Figure  1  shows  a  number  of  data  points.  From  aerodynamic  considerations. 


C[)  is  considered  a  second  order  polynomial  function  of  C|_.  That  is: 

CD  =  ao  +  31  CL  +  32  CL2  (4) 

Using  the  regression  coefficients  a  ,  a2,  and  a2,  the  values  of  C[>  at  the  data 

point  (Ci)  and  the  flight  condition  of  interest  (Ci  )  can  be  found. 

1  2 


% 

"  aO  +  “1  CL!  +  Cl( 

(5) 

cd2 

=  a0  ♦  a,  CLz  *  a2 

(6) 

The 

difference  between  the  values  of 

CD  Is: 

cd2 

CDj  ACD  =  9l  <CL2 

V  * 

a2(CL2  - 

C  2 
Li 

) 

(7) 

cD2 

can  be  redefined  in  terms 

of  CDl 

and  aCq. 

cd2 

=  VaCD  =  CD,  *  a>  (CL2 

•  CLl) 

+  a2  (Ct 

2  _ 

-2 

ck> 

(8) 

Now 

assume 

that 

V 

■  f  <%•  V  V 

(9) 

Using  the  Taylor  series  method  of  error  propagation: 

3Cd2  .  3CD2  3Cd2  2  '-2  do) 

s(cd2)  =  U7T  s(cDi}  1  +  [7r~  sdLl) ]  +trr~s^cL2^]  } 

Dj  LL2 

The  partial  derivatives  in  Equation  (10)  can  easily  be  evaluated  from  the 

definition  of  Cn  in  Equation  (8). 
u2 


They  are: 


3CD2 

"acT" 


(a2  +  2a2  C^) 


(11a) 


(lib) 
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Versus  C,,  Data 


(11c) 


9Cn 

_Z  =  ai  +  2a2  C. 

1_2 

3  C, 

C2 

The  terms  S(Cn  ),  S(C,  ),  and  S(C,  )  in  Equation  (10)  require  interpretation. 

Ui  l i  l2 

They  are  the  precision  indices  of  Cp  and  at  point  1  and  that  of  CL  at  point 
2.  The  CL  at  point  2  is  simply  the  desired  flight  condition.  It  can  be  thought 
of  as  a  mathematical  and  not  a  random  variable  and  can  be  exactly  defined. 
Therefore: 


s(cl2)  =  0 


(12) 


C.  is  a  measured  data  point,  and  because  of  this,  it  is  a  random  variable 

Li 

with  a  mean  and  standard  deviation.  S(C^)  can  be  thought  of  as  the  uncertainty 
associated  with  the  measurement  process.  If  S(C^)  is  constant  or  it  changes 
slowly  in  the  vicinity  of  the  flight  condition  of  interest,  S(C.  )  should  be  a 

Li 

close  approximation  to  the  value  that  S(C. )  would  have  if  C.  were  a  measured 

L  l2 

data  point.  Figure  2  indicates  that  S(CL)  does  not  vary  much  over  the  angle  of 
attack  range  of  interest  (a  =  2  -  6°).  Because  of  this,  S(C^),  where  is 
the  data  point  nearest  the  flight  condition  of  interest,  can  be  used  to  indicate 
the  uncertainty  associated  with  the  measurement  process.  Calculation  of  S(CL) 
is  discussed  in  Section  III. 

Recall  that  the  value  of  Cn  used  in  Equation  (8)  to  define  Cn  was  not  the 

Di  u2 

measured  value  at  C.  ;  it  was  the  value  obtained  from  the  curve  fit  equation. 

Li 

In  other  words: 


\  ■  ao  +  CL,  *  CL* 


(13) 


2 

Therefore,  S ( )  is  actually  the  precision  index  of  ,  S  ,  introduced  by  the 
curve  fit  process.  A  method  for  calculating  S(Cp^)  is  found  in  Section  IV. 

Given  the  values  of  S ( C,  )  and  S(Cn  ),  the  final  calculation  of  the  uncertainty 

L 1  Ui 

in  Cp  at  a  given  flight  condition  can  be  performed  as  follows: 
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U  (CD)  =  tn  S  (14) 

t  is  the  n^  percentile  point  for  the  two  tailed,  student  "t"  distribution. 
Various  percentiles  or  confidence  levels  can  be  chosen.  The  following  discussion 
is  based  on  a  95-percent  confidence  level.  The  ninety-fifth  percentile  point, 
t0.9s»  depends  on  the  number  of  degrees  of  freedom,  which  is  a  measure  of  the 
sample  size  used  to  determine  the  precision  index.  For  S(C^),  t0>95  is  usually 
taken  equal  to  2  since  S(CL )  is  based  on  a  large  sample  size.  This  value  is  not 
appropriate  to  use  for  S(Cg).  The  Aircraft  Compatibility  Branch  uses  a  five-point 
curve  fit  around  the  appropriate  to  find  Cg.  For  this  case,  the  degrees  of 
freedom  are: 

v  =  n  -  k  -  1  =2  (15) 

where: 


n  -  number  of  points  used  =  5 

k  -  number  of  variables  in  regression  equation  =2  (i.e.,  CL#  CL2) 

For  two  degrees  of  freedom: 
to. 95  =  4.303 

The  equation  for  uncertainty  can  be  rewritten  as: 

U(CD)  =  {[ t0  95  S(CD  )]2  +  [t0. 95  (a  +  2a  C.  )  S(CL  )]2}^  (16a) 

'  CD  1  Cl  1  1 

=  {4.3032  S(Cn  )2  +  22[(a  +  2a  C,  )  S(C.  )]2}h  (16b) 

u:  1  2  M  Li 

The  value  U(Cg)  is  the  uncertainty  in  drag  coefficient  at  a  given  flight 
condition  for  a  given  configuration.  More  important  is  the  uncertainty  in  the  drag 
increment  between  the  clean  aircraft  and  the  aircraft  with  the  external  stores. 

The  increment  is  defined  as: 


1C0  ■  CD, 


-  Cr 


Stores  Clean 

From  Reference  1,  the  uncertainty  is: 

UUCJ  =  {U(Cn  )2  +  U(Cn  V) 
u  ^Stores  uClean 


(17) 


(18) 
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SECTION  III 


DATA  POINT  PRECISION  INDEX 

Section  II  described  a  method  to  calculate  drag  data  uncertainties  for 
a  point  on  a  fitted  C^-Cg  curve.  Because  of  the  random  nature  of  the  assumed 
independent  variable,  C^,  the  final  calculation  for  U(Cg)  requires  a  value 
for  the  precision  index  of  CL-  This  can  be  done  in  the  following  manner: 

C^  can  be  defined  as: 

[FN  (cos  a  cosIM  +  sin  a  sinIM)  (19) 

+  FA  (cos  a  sinIM  -  sin  a  cosIM)  ] 

If 

CL  ~  f(xj,  X2,  xn)  (20) 

then  the  Taylor  series  method  of  error  propagation  yields: 

S(C.)  =  {[  S(X!)]2  +  [!5k  S(x2)]2  +  ...  +r  l^L  SU)]*}H  (21) 

dXi  3x2  3xn 

It  is  apparent  that  the  equations  become  more  and  more  complex  as  the 
number  of  independent  parameters  increases.  In  addition,  these  variables 
must  be  independent  or  nearly  so  of  each  other  for  Equation  (21)  to  hold. 
Assumptions  need  to  be  made  on  the  nature  of  the  measurements.  These  are: 

(1)  The  uncertainties  associated  with  some  of  the  independent 
parameters  are  sufficiently  small  that  their  effects  on  the  uncertainty 

of  the  dependent  parameter  are  negligible.  This  is  assumed  to  be  the  case 
for  IM  and  A. 

(2)  The  measured  forces,  FN  and  FA,  are  independent  of  one  another. 

Using  these  two  assumptions,  it  can  be  said: 

CL  =  f  (FN,  a,  FA,  Q)  (22) 

Therefore: 

S(C.)  =  {f^L  S(FN)]2  +  f^L  S(a)]2  +[^L  S(FA)]2  +  {\  S(Q))^ 

3FN  ?o  3  FA  gQ 
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After  evaluating  the  partial  derivatives,  Equation  (23)  may  be  rewritten 
as: 


S(CL)  =  {[(cos  a  cosIM  +  sin  a  sinIM)  ^j^jp'32  +  [Cp  S(<032 

+  [(cos  a  sinIM  -  sin  a  cosIM)  g^-  ]2+  ^CL  S(Q)j2)H  (24) 

For  a  given  test  point,  a,  IM,  Q,  A,  Cg,  and  are  defined.  $(FN),  S(FA), 
S(a) ,  and  S(Q)  are  functions  of  test  conditions  and  instrumentation.  Equations 
for  these  parameters  can  be  developed  in  a  manner  similar  to  that  used  for 
Equation  (24).  $(a)  and  S(Q)  depend  on  the  wind  tunnel  where  the  data  are  taken. 
S(FN)  and  S(FA)  are  dependent  on  the  tunnel  and  model  instrumentation.  Equations 
to  calculate  the  quantities  for  a  specific  case  of  interest  should  be  obtained 
from  the  applicable  test  facility.  The  computer  program  described  in  Section 
V  uses  equations  that  apply  to  the  Aerodynanmic  Wind  Tunnel  (PWT/4T)  at  Arnold 
Engineering  Development  Center  (AEDC). 
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SECTION  IV 


CURVE  FIT  PRECISION  INDEX 

A  method  to  determine  the  precision  index  of  C^  for  use  in  Equation  (16) 
was  just  described  in  Section  III.  The  other  input  required  to  solve  for  U(Cp)  is 
S(CDi).  Recall  that  it  was  defined  as  the  precision  index  of  the  curve  fit  process. 
To  calculate  this  quantity,  certain  assumptions  on  the  nature  of  the  CD  data 
are  made.  These  are: 

(1)  Cp  is  a  random  variable. 

(2)  At  a  given  CL,  possible  values  of  Cp  are  approximately  normally 
distributed. 

(3)  Over  the  CL  range  where  the  curve  fit  is  applied,  S(Cp)  is  nearly 
constant.  This  allows  a  simplified  curve  fit  procedure  to  be  used  (see  pp  106-108, 
Reference  2). 

Since  it  is  possible  to  calculate  S(Cp)  for  a  given  data  point  in  a  manner 
similar  to  that  used  for  S(C^),  the  validity  of  the  third  assumption  can  be 
verified.  Figure  3  shows  a  typical  example.  The  angle  of  attack  region  of 
interest  is  about  2  to  6  degrees.  As  the  figure  indicates,  S(Cp)  does  not  vary 
significantly  over  this  range. 

Recall  from  Section  II,  the  C^-Cp  curve  is  assumed  to  be  of  the  form 

CD  =  a0  +  ai  cl  +  a2  CL2  (25) 

For  a  second  order  curve  fit,  the  normal  equations  to  solve  for  aQ,  aj, 
a2  are,  in  matrix  form: 

[B]  [a]  =  [g]  (26) 

where 
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Index  of  Drag  Coefficient  versus  Angle  of  Attack 


Solution  for  the  regression  coefficients  is  given  by: 

[a]  =  [B]  -1  [g  ] 

Where  [B]"1  is  the  inverse  of  [B],  For  a  five-point  curve  fit,  n  is  equal  to  five. 

In  this  study,  the  value  of  S(Cn)  is  based  on  the  mean  response  vr  r  . 

u  LD|LL,  L2 

The  term  ur  r  may  be  thought  of  as  the  true  value  of  Cn  for  a  given  C,  . 

4)'l'L,  L2  u  L 

In  other  words,  if  many  values  of  Cp  were  measured  at  a  given  Cp,  the  mean  of 

Cn  would  tend  to  ur  . r  r  as  the  sample  size  increased  to  infinity.  Based 
u  ^D1  L,  L2 

on  this,  S(Cp)  is  defined  as: 

S(Cn)  =  s  /[x  ']  [B]-i  [x„  ]  (31) 
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where 


s =  estimated  sample  standard  deviation 

r 

n 


*  [CD.  -  (a0  +  ai  CLi  +  a2  CLi2)]2 


Vs  n 


•FC 


n-k-1 

Cl  23 

lfc 


(32) 


(33) 


(34) 


lFC 

r  2 

FC 

The  FC  subscript  indicates  the  flight  condition  of  interest  where  the  calculation 
is  made.  For  a  more  detailed  discussion  of  the  method  used  to  calculate 
S( Cq) ,  see  a  statistics  text  such  as  Reference  3. 
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SECTION  V 


COMPUTER  PROGRAM  AND  SAMPLE  CALCULATION 

Appendix  A  contains  a  computer  program  developed  in  the  Aircraft  Compatibility 
Branch  to  calculate  uncertainties  in  drag  coefficient  where  CD  is  found  from  a 
fitted  C^-Cp  curve.  Required  inputs  are  explained  in  the  program.  The  uncertainty 
can  be  determined  for  confidence  levels  of  99,  95,  90,  and  80  percent,  i.e., 
x  percent  of  measured  data  points  should  lie  within  the  uncertainty  band  above 
and  below  the  fitted  curve.  The  confidence  level  is  changed  by  usinq  different 
t  values  as  explained  in  Section  II. 

As  noted  earlier,  the  equations  used  to  define  S(CL)  are  applicable  to  the 
four-foot  transonic  wind  tunnel  at  AEDC.  Changes  would  be  required  for  other 
tunnels. 

Assumptions  are  made  to  simplify  the  equations  for  S(CL).  These  include: 

(1)  Only  aircraft  pitch  excursions  are  considered,  i.e.,  6=0. 

(2)  Uncertainty  in  angle  of  attack  is  constant  and  equal  to  0.1  degree. 

(3)  S(a)  equals  one  half  the  uncertainty  in  a  (0.05  degree). 

(4)  A  small  sting  roll  angle  is  assumed  (0.04  degree). 

(5)  Model  weight  does  not  vary  between  configurations. 

(6)  To  a  first  approximation,  model  angle  of  attack  is  equal  to 
sting  pitch  angle. 

(7)  Balance  uncertainties  [S(FNB),  S( FAB) ]  are  functions  of  normal  and 
axial  forces  only  (i.e.,  no  side  load  or  rolling  moment  interactions  are  present). 

(8)  There  is  no  model  roll  angle  relative  to  the  balance. 

(9)  The  precision  indices  of  model  weight  measured  by  axial  and  normal 
force  gages  are  equal  to  the  precision  indices  of  the  balance. 

(10)  Tunnel  total  pressure  is  less  than  1500  lbs/ft2. 

(11)  S(M)  is  constant  and  equal  to  0.002  which  is  40  percent  of  the 


uncertainty. 
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The  assunptions  are  reasonable  for  the  data  of  interest.  Only  pitch  excursions 
are  considered  because  the  performance  problem  concerns  mainly  flight  at  a  constant 
angle  of  attack. 

Appendix  A  contains  sample  calculations  of  U( Cp)  for  a  specific  case.  The 
example  is  for  a  clean  aircraft  at  M  =  0.8  and  CL  =  0.3.  If  the  equation  for 
U(CD)  is  examined  using  the  values  of  S(Cg)  and  S(C^)  from  the  example, 
it  shows  that  most  of  the  uncertainty  in  is  a  result  of  the  curve  fit.  Note 

all  the  confidence  levels  are  shown  in  the  example,  and  U(CD)  decreases  as  the 
confidence  level  is  lowered.  To  be  more  precise,  the  values  S(FN),  S(FA),  S(a), 
and  S(Q)  used  to  find  S(CL)  should  be  redefined  for  each  confidence  level 
instead  of  using  the  95-percent  values.  Since  the  contribution  of  S(C^)  to  U(Cg) 
is  not  large,  this  is  not  significant. 
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SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

A  method  has  been  presented  that  shows  how  drag  coefficient  uncertainties  can 
be  calculated  for  wind  tunnel  data.  Because  of  the  need  to  examine  Cp  as 
a  function  of  C^,  the  uncertainties  in  C^  are  due  both  to  the  uncertainties  in 
CL  and  those  of  the  curve  fit.  While  the  method  of  calculating  the  precision 
index  of  the  curve  fit  is  general,  that  used  to  determine  S(CL)  will  depend  on 
the  wind  tunnel  and  particular  test  instrumentation.  Consultation  with  the 
applicable  test  facility  will  be  necessary  to  work  out  suitable  equations  for 

s(cL). 

The  example  given  in  Appendix  A  indicates  that  the  uncertainty  can  be  fairly 
large  for  a  high  confidence  level.  Reducing  the  confidence  level  results  in  a 
considerably  smaller  uncertainty.  It  would  be  very  useful  if  data  uncertainties 
associated  with  flight  tests  were  given  a  thorough  analysis.  Determination  of 
these  uncertainties  would  increase  the  confidence  in  performance  estimates  based 
on  flight  test  data. 
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COMPUTER  PROGRAM  FOR  UNCERTAINTY  CALCULATION 
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74/74  0s  T  =  3  TPACE 


FTN  4.9*518 


PROGRAM  UNCF9<  IVPJT,0liTPUT,TAPFt;=INP'JT,TAPE6=Cl)TBUn 

THIS  P^GRAM  CALCULATE*  THE  UNCEeTAINTY  IN  DRAG  COEFFICIENT 
*0--  A  SPECIFIC  Flight  CONDITION  OF  INTEREST  ON  A  FITTED 
CL-CP  CURVr.  THF  UNCERTAINTY  IN  CD  IS  HADE  UP  DF  A 
CO“flT NAT  ION  OF  UNCERTAINTY  IN  lift  COEFFICIENT  AND  THF 
UNCERTAINTY  INTRODUCED  bY  CUPVF  FITTING  CL -CO  DATA. 

THE  FOLLOWING  CONSTANTS  ARE  STORED  IN  DLOCFOATA. 

pH 1 1 -A SSUHf 0  STING  POLL  ANGLE  OUPING  A  PITCH  SWEEP  (Of Gl 
SPHI I-STANOAPO  DEVIATION  OF  STING  iOLL  ANGLE  COEGI 
SL-CONSTANT  DEPENDING  ON  PARTICULAR  WINO  TUNNEL 
W-AIPCBAFT  HODFl  WEIGHT  (LOS) 

AG, Ao- CONSTANTS  USEO  TO  OBTAIN  STANDARD  DEVIATION  OF  TUNNEL 
TOTAL  °°ESSUKE .  DEPENDS  ON  TUNNEL  OF  INTERFST. 

Sh-STANOARD  DEVIATION  OF  MACH  NUMftFP 
SAlPHA-STANOARO  deviation  cf  angle  of  attack  <rao> 

STW-ANGLE  nFTwfEN  AIOCPAFT  watfr  linf  and  halance  a X is (OE g > 
A-AIBCBAFT  HCDEL  PFFEBFNCf  AREA  (SO/FT> 

A1 7, A16-C0NSTANTS  USEO  TO  Rf  T  FRM I NF  ETANOARO  PEVIATION  OF 
TUNNf  L  STING  PITCH  ANGLE.  DEPENDS  ON  WINO  TUNNEL. 

Y1.XC.X3  ,Y9-PALANCr  CONSTANTS 

AW1. AK15. SK29-DAL ANCE  CONSTANTS  (CALIFO  K1,K19,<?9» 

SF  AN3,cc  AA2-DALANCF  AXIAL  r0">rr  GAGE  UNCERTAINTIES  WITH 
APPLIED  NORMAL  ANO  AXIAL  LOAD  RESPECTIVELY. 
SFNNE,SFNA?-°ALANCF  NCPHAL  FORCE  GAGE  uncertainties  WITH 
APPLIED  NORMAL  AND  AXIAL  FORCE. 

THE  FOLLOWING  VALUES  APF  CALCULATED  F °CH  The  input 

'■ILPI-STANOARC  DEVIATION  CF  STING  PITCH  tNGLE 
SF A-  "  "  “  AXIAL  FORCE 

SF*_  ••  «•  ••  NOFMAL  FPPCE 

SO-  “  “  DYNAMIC  PRESSURE 

SCL  T  S-  "  “  "  LIFT  COFFPICIFNT 

The  FOLLOWING  INPUTS  APF  *09  THE  TEST  POTNT  CLOSEST  TO  THE 
FLIGHT  CONDITION  OF  INTEREST. 

ALPHA-AIRDFAFT  A  N3l  r  OF  ATTACK 

6-ACH-HACH  N'JMRfo 
PT-TUNNEL  TOTAL  PRESSURE 
D-TUNNEI  OTNAHIC  PRESSURE 
CLTS-LIFT  COEFFICIENT 
C  n  T  S-DR  A  G  COEFFICIENT 

Bpl-VAIUE  OF  LIFT  COFFFICIENT  AT  THE  FLIGHT  CONDITION 
OF  INTEPFST. 

C.Or.F  I  IS  the  confidence  LEVFL  SFLECTFO. 

VALUES  OF  .99.  .°5,  .9,  AND  .  A  ARE  AllOWFO. 

T HE  five  Cl  ANO  CD  VALU*S  ARE  THE  ONES  USED  FOR  THE  CURVE  FIT 

XG.X0P.9.G.  AND  AC  AFf  MATRICES  USED  TO  09TAIN  REGRESSION 

COEFFICIENTS  AND  CALCULATE  THF  S T ANOAfi 0  DEVIATION  OUF  TO  FIT 

AO.Al.A?  A  RF  THE  REGRESSION  COcFFICIENTS  CALCULATED  FO®  THE 
FITTFD  CURVE  WHE  R£  JT  IS  PF  THE  FORM 
CO  =  A’  ♦  A1*CL  *  A2*(CL**£> 

Sc I TE  IS  THf  SOUARE  OF  THE  STANDARD  DEVIATION  DUE  TO  CURVE  FIT 


COMMON/ A /PHI I, SPHI I, SL. W, A5, A 6. SM.S ALPHA, AIM, A, A17, A18 
COMHON/d/X9«£K15,$F AN2.SFAAZ 


0.0  MH  ON/C  /  X  l,XZ,X3,AKl,A<Z9,SFNN2,ScNA2 

cohhon/d/salpi,sfa,sfn,so,sclts 

common/ E/ ALPHA, A  MACH, PT,0, CL TS. COT S.POI 
COMMON/F/CL(5l,CD(5).XO(3»,XOP<3»,DC3,i», 


G  (  3  )  .  AC  ( 3 ) 
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frocram  wmcf* 


TMfJ 4  ORT*0  RRACE 


FT*  4.* *510 


CONN  Off/  IT/  AO  •  *1  •  A775TI TT - - - 

C0NN0N/M/A0<U,  AL « 3 1  . AMI  3) .Z( 3) . X(l> 

18  RE  AO  <»  .1 ) AL^HAt ANACH*PI»Q»CLTS,COTStPOI - 

1  FORMAT (TF10.  5 ) 

IF (EOF (5)  1 99*  3 

1  READ'S, 81  -eMfl  —  -  —  - 

a  FORMAT {F10  >51 

mm\  miiHHHii - 

2  FORMAT (5F10. 51 

CALL  SOCL  - - 

CALL  SOFIT 

WRIT  E(6, 111  AHACH 

11  FORMAT  <lHl,t  JHMAet^  NUMBER  *,Flfl.81 - -  - 

NRITE (6,12) 

14  atHu^awHSiMF1  »«* - 

14  FORMAT I/5X.26MVALUES  OF  DRAG  COEFFICIENT/1HO ,5F12 .4) 


WRITE (6, 15  >  POI 
FORMAT  t/F/KX. SOHttFT 
WRITE (6.  71  AO . A1 , A2 


IK  9 - 


COEFFICIENT  OF 


115 _ 


1  &f if E  (6*  4*’  «?^lzJSlE 

_  __4 _ pqpmaT t/'lFiZGHSTAMDARO  OEVIATIOH  |M  CL  *iF1TiM 

SFIT*SQRT(SFIT2» 

_ MRITET6.5>bSFIT  DEVIATION 

9  CONF»?ONFI*100?ST*,<OAtO^  !OM  IK  fM-r.TtMrt 

1 6  *  x ^?HC3hFI6EMCE  INTERVAL  *,r5.8.ZHTn 

_ lfIgS®W8;i;8UJjO|j _ 

IF(CONFI.EQ.0.9>  GO  TO  tiJ 
IF (COMF1 .EQ.0.8)  30  TO  103 

-18(T--T1T*Lt»r.325  - 

TVAL 2*2.58 

TTT-fvA^ilf^J - 

.  _  -ins-s« _ 

182  TVAL 1*2. 92 

TVAL 2=1, 65  _ 

193  iw  1*1.886 

7tt  ~  VAR1  -  Hv?L  I  *  *2  .T  ▼sFl  rs - 

VAR2«(T!fAL2**2.»*((AlF2.*A2*CLTSI**2.»*(SCLTS**2.» 

- UCO*  3QRT1V-AR1  »VAR2  ♦  — - - 

WRITEI6.6I  UCO 

6  FORMAT (///IX. 3 TH THE  UNCERTAINTY  IN  DRAG  COEFFICIENT 

- - 

_ £Nfl-_  _ _ _ 


.F12.4I 
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iif  gam 

3UUKUU1  lrlt  JUt/t 


7M/T+ -  0PT*6  TRACE 


— j_ 

5 

10 

15 

20 

2? 


X 

-i- 

c 

■4 


SUBROUTINE  SOCL 

^HIS_ SUBROUTINE  CALCULATES  -SCtTS 


LIFT  COEFFICIENT  FOR  THE 
CONDITION  OF  I NTEREST . 


CnS2°H/*/?q1lKlgHsr)|N^^M?*&’SMtS*LPM*’*IM’**  Al7’*18 
COMNON/C/X1* X2, X 3, A K1 ♦AK29 ,SF  NN2 iSFN»Y 
C0NH0N/0/S«LPI,SFA,SFN.SQ,SCLTS 

- C8HN0N/E/ALPHA. AMACHvPTv8*CLt5»C0T S»PQI  -  -  —  — 

SALPI*A17*A18*AL PH* 

_ L?F  AN2* *1,  *SEiAf*i  t  >  I _ _ 


SFAG2=ISL*X9*AK15>**2.*SFAB**2. 

SFAST2* lM*COSO (ALPHA! *SALPH  * *2. * (SINO  C ALPHA! *SHA ) • *2. 
SFA*  <SFAG2*SFA$T2»**.5 
SALPHA*. 000873 

9PNfr*.5»  HSF  NN2**E . ♦SFNAE**E.+***Et 


SFNG2=USL/Xl»**2.)*t  (X2*AK1)  **2 .»<  X3*AK29I  **2.» 

_1 _ *SFN0**2 

SVOUSFNB 


3S 


__«l» _ 


l/ARl  N*  (N*SIND(AL  PH  A 1  •COSO  TPHII »*SALP| »**2. 

.  _  VAR2N=t COSOt ALPHA!  ’COSCIPMII1* - 

VAR3N*(M*C0S0(AlPHAI*SINDCPHIII*SPHII»**2. 

SFNST2=VAR1N*VAR2N*VAR3N 

— SfN*TSFNS2*SP-NST2»**iS - - - 

WAR1P* ( ( f2.*0/AHACHI*ll.-1.2* ( AMACH**2.  II  »  / 
1  (1.*.2*<AHACH**2.  )M*SN 

SPT*1T*A6»PT - J - *“ - 

l/AR2P=(Q/PT»*SPT 
—  SO* tVARLP**2.*WAR2P**2.J**» 


WAR1T* ((COSO (ALPHA! *C0SDfAIN>  *SI NO I ALPHA I *S I  NO I A INI  I 
1  *SFN/IQ*An**2. 

UARST* (COT 5*SAL PHOT  - — - — - 

VAR3T*lTC0SfMALPHA>*SIN0UIM»-SINOUlPHAt*COSOUIH»» 

1  •SF»/(0*A» 1**3 


SCyJS=SQRTmRlT*«R2T*VAR3T*VARAT» 

EN0 
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t 


SUBROUTINE  SOFIT 


74/74  OPTrfl  TRACE 


FTN  4.0P510 


10 


15 


20 


25 


50 


»- 


Jt9 


45 


50 


59 


60 


C 

c 

c 

c 

e 


SUBROUTINE  SOFIT 

THIS  SUBROUTINE  CAtCULATfS  -THE  REGRESSION  COEFFtOifNTS  TO  FIT 
THE  CL-C0  OATA  AND  DETERMINES  THE  STANOARO  DEVIATION  OF 
.THE  CURVE  FJT,  _  _ _ 

COMMON/E/ALPHA,AMACH,PT,Q,CLTS,CDTS.POI 

COHNON/F /CL  ( 5) .CO  15 1 . XO  <31  ♦X»PESA-,ali.3U«Sl,AC433 -  — 


_.c _ 


COHNON/G/AO,  A1.A2.SFIT2 
COMMON/H/AO(lt  .ALUI.AMUI  .ZU).X(l) 
S  UH1 1*6  • 

SUM? 1=0. 


10 


29" 


O 

C 

J 


SifMll=0. 

SUH4 1=0. 

00  10  1=1,5 
SUM11=SUH11*CL  Ill 
SUM21=SUH21*CL(II»*2. 
SUH31=SUiSl*Ct II»**3. 
SUM41=SUH41*CL CII**4. 

-mm.-  —  — 

0 ( 1. 2) =SUM11 
B  11.3) &SUH21 
B  ( 2, 1I=SUM11 
B (2. 21 =SUM21 
Bt2.5»=SUM31 
B(  3, II  =SUM21 

IfiliUiifflft 

i!KS§:8: 

SUH3  2=0 . 

DO  20  1*1,5 
SUH12*SUH12»C0«II 
SUH22=SUM22»CL (I l*CDII) 

mm**'-'*-™-*  *-?JE0  >L> 

Gill =SUH12 

G (?) *SUH22  -  - 

GUI  =SUM32 

CALL  MINV(B,3,AD, AL.AMI 


CALL  TO  SYSTEM  ROUTINE  TO  INVERT  MATRIX  8 


call 


MULTIPLY  INVERSE  OF  fl  BY  G- SO -OBTAIN  MATRIX  AC 
A0=A  C(l) 

A1=AC(2) 

A2  =  ACtSI-  - - 

SUH3=0. 


Sum! ^CT)  TTT -A 0-Al*CL II J- A2*  fCL 111**2.1 1 

1  **2. 

CONTINUE  - - - - 

SSO*SUM3/2. 


X0PUI*1. 
RJft?T*fOI. 


W 


TO 


40 


C 

£ 


XfiPt 31 f POI **2 . 

CONf INU 
■CAtr 


RtMKOP.B.Z.1.3,31 


JWLTIPLY  XOP  BY  INVERSE  OF  B  TO  OBTAIN  Z 
CALL  GMPRO(Z,XO«X,1,3.1I 


— 75 


$ — ■RDtrpLrZ  ftY^  jffl  TO  OffTAIN  Y7  X  IS  A  OWE  COLUMN  OWERDM^ 


23 


V 


SUBROUTINE  SO* IT -  7*r/M  OBT~fr  TR*CE  - 


trvrarro>r  doe  to  curve- tit. 

XlaXtll  -  - - 

SFIT2=SS0*X1 

RETURN 

END  -  - .  - 


— BtOCK  OAT*  BtKOAT. - 74 /T4  — 0*T-=tt  f**et - PTH  4*8+818 


' 5  ■  1 

BCOuRQATA 

ESKS'/^^lil^sfU 

jA6,5N,S ALPHA, AIH.A.A17, A1B 

COHHOM/C/Xl. X2.X3. AK1.AK29,! 
DATA  PH{I,SPHXItSL,MtA5,A6,: 

SFNN2.SFNA2 

SM,$ALPHA,AI«,A, A17.A18/ 

•  in  l.otoo,  «  «i?f*uuuH««DUCM  vouar  u« - 

• 75t« 0 04 t» 00 0 14/ 

DAT  A  X9«  AK15  •SFAN2,  SF  AA£/l_*_i«0  359*«4i»45/  _ 

DATA  Xt,XZ,Xr,A<l,AK29,SFNNZ,SFNAZ/3.91Z,-l.,l., 

10  1  .303, .3994. .89, .76/ 

-  . . -6MO . -  . 

MINO  TUNNEL  DATA  POINTS 

VALUES  OF  LIFT  COEFFICIENT- 

- 231 —  -  .3324 .V325  - »«Wi» 

_ VALUES  OF  OPAG  COEFFICIENT _ _ 

.0152  .0163  .0272  .046*  .0603 

—  -t  I  FT— COEFFICIENT  ©F  I-HTE*EST  ■  *300©  - 


CURVE  FIT  COEFFICIENTS 


AO  = 
A1  * 
A  2  » 


•01664 

••04606 

.22537 


ST  ANOARD  OFVTrriON  IN  "Cl  “* 


STANDARO  DEVIATION  IN  FIT  = 


TPT!3T' 

.  00 14 


confidence  interval  =  P9.  y. 


-IH £  UNJCER.TAWTY  IN  DRAG  COEFFICIENT 


.£139 _ 


nAC“  mJnotrt  -  • 


MIND  TUNNEL  DATA  POINTS 

VALUES  OF  LIFT  COEF FIC IFNT 

- -y -  vtT^t  *3324-  - —  *4*25 

_ VALUES  OF  DRAG  COEFFICIENT _ 

•0152  .0169  *0272  .0466 

- LIFT  COEFFICIENT  OF  INTEREST  = - .3003  - 


CURVE  FIT  COEFFICIENTS 


.01684 

>.04606 

v225*F 


“  ^TTWDATrtrorVTA^  CTT~=  TOOTS 

STANDARD  DEVIATION  IN  FIT  =  .0014 


CONFIDENCE  INTERVAL  *  95.  7. 


THE  UNCERTAINTY  IN  DR  AO  COEFF ICIENT  = 
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HAC*  NUMBER 


NINO  TUNNEL  DATA  POINTS 

VALUES" OF  LIFT ■COEFFICIENT 

- *0231 - rt-W - - *4925 

VALUES  OF  DRAG  COEFFICIENT _ 

•0192  *0169  .0272  .0468 


H.IFT  -COEFFICIENT  OF -INTEREST-  =  -  .-3989- 


CURVE  FIT  COEFFICIENTS 


AO  * 
At  = 
A2~* 


.01684 
-• 04606 
.-22547- 


WOIRO  DEVIATION  "ITT  TL  = 
ST AMD A 90  DEVIATION  Tn  FIT  = 


CONFIDENCE  INTERVAL  *  90.  Z 


.inr 

•  0014 


6474 


0803 


THE  UNCERTAINTY  !!*  QRAG— CQEFF  IC-IENT-* _  .C041 
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V 


- NftCH  NUN8SR  - .700  - 


KIND  TUNNEL  DATA  POINTS 

VALUES  ”oF~LiFT "COEFFICIENT 

- - rl776 - .JJ24 

_ VALUES  OF  DRAG  COEFFICIENT _ 

.0152  .0169  .0272 


— — LIFT  COEFFICIENT  OF  I  NTFRES7  -  —  .3008 - 


CURVE  FIT  COEFFICIENTS 

AO  *  .01684 

- . 

STANDARD  WVTATIffN  TflTCt  * - .0033 

STANDARD  DEVIATION  Tn  FIT  =  .0 01 4 

CONFIDENCE  INTERVAL  =  80.  X 


- TME  UMCERTA3NTV  IN  CRAG  COEFFICIENT  =  - .0027 
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